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ABSTRACT 

We analyze the temporal evolution of the three-dimensional (3D) magnetic 
structure of the flaring active region (AR) NOAA 10930 by using the nonlin- 
ear force-free fields extrapolated from the photospheric vector magnetic fields 
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observed by the Solar Optical Telescope on board Hinode. This AR consisted 
mainly of two types of twisted magnetic field lines: One has a strong negative 
(left-handed) twist due to the counterclockwise motion of the positive sunspot 
and is rooted in the regions of both polarities in the sunspot at a considerable 
distance from the polarity inversion line (PIL). In the flare phase, dramatic mag- 
netic reconnection occurs in those negatively twisted lines in which the absolute 
value of the twist is greater than a half-turn. The other type consists of both 
positively and negatively twisted field lines formed relatively close to the PIL 
between two sunspots. A strong Call image began to brighten in this region 
of mixed polarity, in which the positively twisted field lines were found to be 
injected within one day across the pre-existing negatively twisted region, along 
which strong currents were embedded. Consequently, the central region near 
the PIL distributed with a mix of differently twisted field lines and the strong 
currents may play a prominent role in flare onset. 



Introduction 



Solar flares and coronal mass ejections (CMEs) are eruptive liberation of the accumu- 
lated free magnetic energy in the solar corona and are considered the biggest explosions in 
our heliosphere. These phenomena affect geospace in the form of electromagnetic distur- 
bances called geomagnetic storms. Therefore, it is an important issue for the space weather 
forecast to have the better understanding of the physical reason responsible for triggering 
these phenomena. 

Many models based on the magnetohydrodynamic (MHD) approach have been proposed 
and have indicated that twist and shear of the magnetic field lines in particular are responsible 
for these processes. The twist number of the magnetic field is crucial for analyzing the 
stability of the magnetic flux tubes. For instance, in a periodic system such as a cylinder, 
a twist with more than one turn could desta bilize a flux tube; this is widely known as the 



Kruskal-Shafranov limit f lKruskal fc Kulsrud 



19581). The dyna mics of cylindr i cal flu x tubes 



in the solar corona have been investigated by lForbesI (jl990l ) and llsenberg et al.l ( 1l993l ) in two- 
dir nensional (2D) space, a nd their work was later extended to three-dimensional (3D) space 
by llnoue &: Kusand (120061 ). In an anchored flux t ube, such as a coro nal loop, a stronger twist 
is required to destabilize the ideal MHD modes ( iTorok et al.ll2004j ). Therefore, an accurate 
quantification of the magnetic twist is necessary in order to clarify the role of the ideal MHD 
instabilities in triggering solar flares. However, there is still no consensus on the question of 
how strong the twist and shear of the magnetic field hues must be to trigger a flare. 



- 3 - 



Active region (AR) NOAA 10930 produced an X3.4-class flare at 02:10 UT on Decem- 
ber 13, 2006, and also gene rated a coronal r nass ejection (CME) that caused electromagnetic 
disturbances in geospace ( iLiu et all l2008l : iKataoka et al.ll2009l ). Flare-associated features 
(e.g., the flare r ibbon. X-ray sigmo id, and cusp loop structure) in this region were observed 
well by Hinode ( iKosugi et al.ll2007l ). Furthermore, Hinode successfully conducted continuous 
observations of the photospheric magnetic fleld corresponding to this AR at an unprecedent- 
edly high resolution before and after the X3.4-class flare occurred. Therefore, this is an ideal 
object to study in order to understand the characteristics of ARs producing X-class flares. 

The positive sunspot in AR 10930 was associated with strong sheared and twisted 
counterclockwise motion, whereas the negative sunspot was almost stationary compared 
to the positive one. As a result of the counterclockwise motion of the positive s unspot, a 



negati ve (left-handed) twist was injected into the overlying coronal magnetic loops. ISu et al. 



(120071 ) reported a strong shear fleld on the basis of X-ray observations in the solar corona. 
The development of the magnetic energy and its injection into the corona in the form of 
magnetic helicity becaus e of the strong shearing a n d twisting moti ons o f the sunsp o t wer e 
also reported recently by lMagara fc Tsunetal (120081 ) . ISu et al.l (120091 ) . and lPark et all ( 120101 ). 



The apparent motions of the po sitive sunspot asso ciated with this region were also 



reported in some previous studies; e.q.. lZhang et al.l ( 120071 ) investigated the sunspot rotation 



using white-light images from the Transition Region And Coronal Explorer and estimated 
its rotation for thre e days beginn i ng on December 11, 2006, flnding that it was about 240°. 
On the other hand, iMin fc Chad (120091 ) estimated that it was about 540° by using G-ban d 
images taken by the Solar Optical Telescope (SOT) on board Hinode (ITsuneta et al.ll2008l ). 
Furthermore, the local area surrounding the polarity inversion line (PIL) contains small- 
scale complex magnetic structures of mixed polarities due to the twisting motion of the 
positive sunspot. Some authors al so suggested the formation of such structures because 



of flux ernergen ce around the PIL ( IZhang et al.l 120071 : iKubo et al.l 120071 : IWang et al. 



Lim et al 



2010 



Ravindra et al 



) and indicated that it plays a key role in triggering flares 



Park et al 



20081 : 



20101 : 



201ll ). It is obvious from the above information that the apparent rotation 



of the sunspot might be rather ambiguous for quantitative measurement of the magnetic 
fleld twist. Therefore, it is difficult to obtain appropriate information on the fleld line twist 
in 3D space. 

Consequently, nonlinear force-free fleld ( NLFFF) extrapolation becomes a powerful tool 
for understanding the 3D magnetic structure. ISchrijver et al.l (120081 ) applied NLFFF extrap- 
olation to AR 10930 and identifled a strong electric current regio n above the n eutra l line 



before the flare, most of which disappeared as the flare proceeded. Ilnoue et al.l ( 20081) also 



indicated the possibility of fleld line relaxation during this flare. Furthermore, 



Guo et al. 
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( 120081 ) suggested that a magnetic dip might be form ed in the AR and suggested that it could 
sustain a filament above the magnetic neutral line. iJing et al.l (|2008[ ) analyzed the altitude 
variation in the magnetic structure in the pre- and post-flare phases and found that the 
energy release process proceeded in some height range from ~ 8 Mm to ~ 70 Mm, whereas 
the non-potentiality of the magnetic field increased after the flare below ~ 8 Mm. However, 
there is still no commonly accepted explanation for the entire flare dynamics {i.e., from 
energy accumulation to relaxation) based on the quantitative twist or topology of the field 
lines. 



Recently, llnoue et al.l (1201 if ) also developed an NLFFF extrapolation procedure based 
on the MHD relaxation method and applied it to AR 10930. They confirmed its reliability 
by comparing the location of the footpoints of the sheared field lines across the PIL with that 
of the Call illumination obtained by SOT / Hinode. Their results show that the footpoints 
before the flare correspond well to the location of the Call illumination in the central area 
of the entire domain. They also introduced the magnetic twist, which represents the degree 
of twist in the magnetic field lines, and clarified that the many strongly twisted lines in this 
AR have less than a one-turn twist, which indicates robustness against kink instability. The 
fraction of the magnetic flux corresponding to the strongly twisted lines having more than 
a one-turn twist was found to be negligibly small compared to those of less twisted (^0.5- 
to ~ 1.0-turn twist) field lines. On the other hand, in a later study, llnoue et al.l (120 12[ ) 
investigated the 3D magnetic structure making up the sigmoid and clarified the relationship 
between the X-ray intensity obtained from the X-Ray Telescope (XRT) on Hinode and the 
magnetic twist or field-aligned current obtained from the NLFFF. They found that strong 
X-ray intensity is closely related to the field-aligned current flowing in the chromosphere 
rather than the twist values of the strongly twisted field lines. They further indicated that 
the field line patterns generated by the NLFFF a re quite similar to the profiles obtained 
in the flux-emergence simulation of iMagaral (120041 ). In addition, other topics, such as the 
characteristics of the 3D magnetic field in this AR, the formation process of flare-producing 
ARs, and a quantitative interpretation of the flare dynamics in terms of the magnetic twist, 
were not addressed in our previous studies. 

In this paper, we also analyze the magnetic twist of AR 10930 by applying NLFFF 
extrapolation to time series vector magnetogram data obtained from Hinode in order to 
understand these problems. The rest of the paper is structured as follows. The numerical 
model and data set are described in Section 2. The results of our analysis of the 3D structure, 
magnetic twist, and topology of AR 10930 are presented in Section 3. Their implications for 
the flare onset mechanism are discussed in Section 4, and the key conclusions are summarized 
in Section 5. 
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Numerical Method 



2.1. Extrapolation of the NLFFF 



We adopt th e sa me force-free extr apolation method as in our previous studies {e.g., 
Inoue et al.ll201ll and llnoue et al.ll2012l ). which employed the MHD relaxation method. A 
3D field ideally has to be extrapolated so that it satisfies a force-free condition, such as 

VxB = aB. (1) 

Unfortunately, the photospheric field o btained frorn an obs ervat ion is not necessary to satisfy 
the force-free condition. Nevertheless, llnoue et al.l ( l201ll ) and llnoue et al.l (120121 ) carefully 
evaluated the reliability of an extrapolated field by comparing it with Call and X-ray images 
from the SOT and XR T on board Hinq d e. They found that their resul t s are q uite consistent 



with the observations. ISu et al.l (120091 ) . ISavcheva &: van Ballegooijenl (120091 ) . and iHe et al. 



(120111 ) also compared their force-free field with multi-wave observations. The magnetic field 
lines obtained from their results were also found to be quite reasonable for recapturing the 
X-ray or EUV images obtained by Hinode; hence, these results support the suggestion that 
the NLFFF is a solid tool for describing 3D magnetic structure. 

In this study, the solved equations are governed by the MHD-like equations for a low-/3 
plasma. 



dv 1 

— = -(v ■ V)v + - J X B + uV^v, 

at p 



dB 

~dt 



V X X B - 7] J) - 



(2) 
(3) 



V X 



(4) 



(5) 



where the last equation was originally introduced by iDedner et al.l ( l2002l ) for calculating an 
MHD solution that also avoids deviation from the solenoidal condition V • S = 0. Here, 
B is the magnetic fiux density, v is the velocity, J is the electric current density, p is 
the pseudo-density (which is assumed to be proportional to |-B|), and is the convenient 
potential. The length, magnetic field, velocity, time, and electric current density are nor- 
malized by Lo = 5.325 x 10^ (cm), = 3957 (G), Va = 5o/(/ioPo)^/^ = Lq/Va, and 
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Jo = Bo/fioLo, respectively. The non-dimensional viscosity u is fixed at (1.0 x 10 and 



the no n-di mensional res i stivity t] is assumed to have the functional form given in llnoue et al. 
J201lh and llnoue et"aD tom . 



V = Vo + Vi 



\J X B\ 


\v 


2 




B\ 



(6) 



where r/o = 5 x 10 ^, and ?7i = 1.0xlO ^, in non-dimensional units. The other parameters, 
c1 and Cp, are fixed at 0.04 and 0.1, respectively. 

The magnetogram set on the bottom boundary is a hybrid map from SOT /Hinode and 
the Michelson Doppl er Imager (MDI) on the Solar and Heliospheric Observatory [SOHO) 
(IScherrer et al.lll995[). T o obtain an extensive overview of the numerical box (see Figure 



Ha) i 



m 



Inoueetal 



20121 ). the magnetogram from the spectropolarimeter (SP) on Hinode 
is located at the center of the MDI/ SOHO data. The area outside of the region covered 
by the SP magnetogram is maintained by the longitudinal field of MDI/ SOHO, whereas its 
tangential components are fixed by the potential field derived on the basis of a synoptic chart 
of the line-of-sight component of the magnetic field observed by MDI/ SOHO. The lateral 
and upper boundary conditions are also fixed by the potential field components. 

An initial condition is given by the potential field calculated from the normal component 
of the magnetic field on all the boundaries after satisfying J BndS = in the entire domain, 
where dS represents the surface element on all the boundaries, and the subscript n represents 
the component normal to the sur faces of the boundaries. The 3D configuration is shown in 
Figure |2](b) of llnoue et al.l (120121 ). The velocity field v set to zero on all the boundaries. 
A Neumann-type boundary condition ((9„0 = 0) is applied to the potential at all the 
boundaries, where dn represents the derivative along the normal direction on the surface. 

The numerical scheme for this calculation is given by the Runge-Kutta-Gill scheme with 
fourth-order accuracy for the temporal integral and the central finite difference with second- 
order accuracy for the spatial derivative. The simulation domain is a rectangular box span- 
ning (0,0,0) < {x,y,z) < (4.0Lo,4.4Lo,2.2Lo), which corresponds to (295.2", 324.6", 162.3") 
in view angle. The domain is uniformly divided by a 128 x 128 x 64 grid. The vector- field 
magnetogram located on the bottom boundary (125 x 64 grid) was formed by binning from 
the original magnetogram, which had a 1000 x 512 grid. 



2.2. Observations 



We extrapolated the 3D force-free field of AR 10930 by using the vector magnetograms of 
an observation provided by the SP on SOT / Hinode. We used five magnetograms, which were 
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observed at 17:00 UT on December 11; at 03:50 UT, 17:40 UT, and 20:30 UT on December 
12; and at 04:30 UT on December 13, 2006. The first four data sets were taken before the 
onset of the X3.4-class flare, which occurred at 02:10 UT on December 13, 2006, whereas the 
last one was observed just after it. These data were obtained by Milne-Eddington inversion 
of the Fel lines at 630.15 nr n and 630.25 nm. The minimum energy method w as apphed 
to solve the 180° ambiguity jMetcalj \19M . iMetcalf et alll2006[ Iheka et allboooh . We also 
compared the Call H 3968 A images with the 3D NLFFF. These images were provided by 
the Broadband Filter Imager (BFI) on SOT /Hinode, which has a time cadence of 2 min. 
The field of view is 223.15" x 111.58" , with a pixel size of about 0.109". The Call line is 
sensitive to temperature at lO^K, which corresponds to the lower chromosphere, and reacts 
strongly to chromospheric heating. In this study, the Call images were employed from just 
after the flare onset to the growth phase of the two-ribbon structure, which corresponds to 
the time interval from 02:14 UT to 02:40 UT on December 13. 



2.3. Analysis of the magnetic twist 



We focus on the magnetic twist, which is defined as the turnover number of magnetic 
field lines in magnetic flux tubes and plays a key role in judging the stability or instability of 
a magnetic configuration, as mentioned in section 1. Although the detailed formulation and 
its implementation in terms of the results obtained from t he NLFFF on AR 109 30 using this 



magn etic twist are illustrated in our previous works, (e.q.. llnoue et al.ll201ll and llnoue et al. 



2OI2I ). we briefly describe this key issue again here to refresh readers memories. 



The magnetic twist (T„) is related to the magnetic helicity; e.g., the hel icity of a closed 
flux tube anchored on the solar surface is described by the following equation ( iBerger fc Field 
19841 : iMoffatt fc Riccalll992l : iBerger fc Prioi] 120061 ) . 



H={T^ + Wr)^' 



(7) 



where H, $, and Wr are the magnetic helicity, magnetic flux of a cross section of the flux 
tube, and magnetic writhe corresponding to the helical structure of an axis of the field line, 
respectively. T„ indicates how much of the magnetic helicity is generated by the currents 
parallel to the flux tube axes Jm; thus, T„ is described by 



T 

■J- n 



ds 



-ds 



J\\ 



-ds. 



where 1 1 indicates the component parallel to the field line , and the line iiitegral f ds is taken 



along the central magnetic field line of the flux tube (IBerger fc Priori l2006l : iTorok et al. 
2010h . 
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In this study, because we assume the force-free state and calculate T„ for each field line, 
the twist value can be replaced by 

If 1 

Tn = — / adl = — aL, (9) 
47r J An 

where a is the force- free a given in equation ([T]), the line integral J dl is taken along the 
magnetic field line, and L represents the field line length. 

We introduce the average force- free a, which is defined as 

a = ^(a+ + «-), (10) 

where a+ and a_ indicate the force-free a on the footpoints of each field line on the positive 
and negative polarities, respectively. This formulation is adopted to minimize the deviation 
of the force- free condition in the photosphere. Consequently, the twist is formulated in this 
study as 

Tn = ^aL. (11) 

Before calculating the average force-free a (a), we first focus on those specific magnetic fields 
for which the normal component exceeds 30 G in order to avoid numerical errors; further, 
a = is assumed for the rest of the area. Next, we average a over 3x3 cells, and finally we 
obtain the average force-free a {a) after calculating and 

Note that we should use this formulation of the twist with caution, as given in equations 
(IH]) or (ini). Because its value also depends on the field line length, one would obtain very 
large values corresponding to the very long field lines that have one footpoint outside the 
AR. However, in this study we do not consider such long field lines and essentially restrict 
our analysis to closed field lines, using the original definition in which the footpoints of 
the magnetic field lines are anchored somewhere inside the AR, and the magnetic twist is 
considered as being effectively propagated along the field lines through the strong shearing 
and twisting motions of the sunspot. 



3. Results 

3.1. 3D structure of the Solar Active Region 10930 

First, we try to understand the characteristics of the 3D magnetic structure 6 h before 
the flare associated with AR 10930 using twist analysis. The choice of this time is obvious 
because the vector field in its final pre-fiare phase was observed around this time. The details 
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of the force-freeness of the extrapolated fields were given in our previous work, llnoue et al. 



( 120121 ) ■ Figure [Uj^a) shows the distribution of the normal component of the magnetic field 
with the selected contours in white solid lines corresponding to = 790 G and —790 G, 
whereas the blue dotted lines correspond to the PIL. The black solid square box indicates 
the region surrounding the PIL between the positive and negative polarities of the sunspot. 

Figure [U^b) represents the twist distribution of the field lines of the NLFFF mapped 
on the photosphere in the same field of view as in Figure [T]^a), where the white solid and 
blue dotted lines have the same meaning. Note that these twisted values are focused on 
closed field lines only, as mentioned in subsection 2.3. The positive and negative twist 
distributions describe right- and left-handed field line twists, respectively. A strong negative 
twist (|T„| > 0.5) is clearly distributed on both sunspots; the strongest twisted region (more 
than one turn) also appears here. Hereafter, "strong twisted lines" usually refers to field 
lines having more than a half-turn twist (|T„| > 0.5). From this result, it is evident that the 
energy accumulated field lines are connected mainly between these strong twisted regions, 
which are at a considerable distance from the PIL. On the other hand, the region near the 
PIL surrounded by the black solid square in Figure [D^a) is of mixed polarity, containing both 
positively and negatively twisted field lines. However, the magnitudes of the twist values 
are much smaller (|T„| < 0.25) than those of the outer field lines characterized by strong 
negative twist that occupy the main twisted region of both polarities 

Figure [T|^c) shows the field line length (L) mapped on the photosphere in the same 
field of view as Figure (U^a), where the regions surrounded by red lines correspond to a twist 
value of \Tn\ > 0.5. The region marked by diagonal lines is occupied by open field lines 
that are connected outside the field of view. The mixed region of positively and negatively 
twisted field lines near the PIL enclosed by the solid square in Figured]^ a) consists of shorter 
magnetic field lines compared to the strongly twisted region surrounded by the red contours. 
From these results, AR 10930 can be roughly divided into two regions. One is occupied 
by the strongly negatively twisted field lines of single magnetic helicity and is located at 
a considerable distance from the PIL. The other corresponds to a mixed region of both 
positively and negatively twisted field lines near the PIL between the positive and negative 
polarities of the sunspot. 

A 3D representation of the magnetic field lines is shown in Figure [U^d), which covers 
the same field of view as Figure [TJ^a). The blue, orange, and green field lines represent 
selected 3D magnetic field lines having twist values of |T„| < 0.5, 0.5 < |T„| < 1.0, and 
1.0 < \Tn\ < 1.5, respectively. The orange and green lines are clearly longer than the blue 
one. On the other hand, the magnetic shear corresponding to the blue lines is stronger than 
that of the other field lines despite its weaker twist strength. This will be discussed in more 
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detail in a later section. 



3.2. An Evolution of the Magnetic Twist in the Solar Active Region 10930 

3.2.1. An Evolution of the Strong Twisted field Lines Forming at a Considerable Distance 

from the PIL 

In the previous section, we investigated the 3D magnetic structure of AR NOAA 10930 
6 h before the flare. Here, we investigate the temporal evolution of the magnetic twist 
associated with this AR using the NLFFF derived from the time series of the vector flelds 
taken before and after the flare. In the upper panels of Figure [2]^a—e), the temporal evolution 
of twist having three different values, |T„| = 0.5, |T„| = 1.0, and |T„| = 1.5, is plotted in 
red, green, and blue contours, respectively, over the normal component of the magnetic fleld, 
represented in gray scale. The fleld of view is the same as in Figure [1] In the lower panels 
(a'— e'), the orange and green lines represent selected 3D magnetic fleld lines having twist 
strengths of 0.5 < |T„| < 1.0 and 1.0 < |T„| < 1.5, respectively. 

At 17:00 UT on December 11, the strong twist > 0.5) was distributed on the 

positive sunspot and the southwest edge of the negative sunspot. It was built up by the 
counterclockwise motion of the positive sunspot, as mentioned in section 1. The strongly 
twisted regions are connected with orange fleld lines, as shown in the Figure |2]^a'). At 03:50 
UT on December 12, the twist distribution grew wider, and in some regions, even stronger 
twist (more than one turn) can be seen. A comparison of the twist analysis at three different 
times before the flare, 03:50 UT, 17:40 UT, and 20:30 UT on December 12, reveals that the 
regions having a twist of more than one turn were spreading more at both polarities, and 
some magnetic lines with strong twist appear in the expansion stage owing to the continuous 
shearing and twisting motions of the positive sunspot. Figure |2]^d), obtained 6 h before the 
flare, is the flnal picture of the magnetic conflguration in the pre-flare phase. Although the 
strongly twisted lines (|T„| > 1.0) are widely distributed compared to the previous times, 
they are localized on the edge of the sunspot, where the magnetic fleld strength is weaker 
than at the central part of the sunspot. The post-flare image. Figure [2I^e) , corresponds to 2 
h after flare on set. The twisted fleld lines decrease, which is also reported for another AR in 



a recent paper ISun et al.l (120121 ). in the central part of the positive sunspot, which consisted 
mainly of fleld lines of twist |T„| < 1.0 before the flare. However, the part of the strong twist 
conflned at the edge of the positive sunspot remained even after the flare. 

From these analyses, we infer that although the positive sunspot is capable of generating 
strongly twisted lines of more than half a turn in the energy accumulation process because 
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of its strong twisting and shearing motions, the strongest twist (in particular, |T„| > 1.0) 
does not seem to contribute efficiently to the main energy release. To further investigate 
this point, let us examine the distribution of the current density along the twisted field lines. 
An estimation of the strong current density, which is related to the measurement of a large 
free energy through the twisted field lines, would provide a way to clarify whether the most 
strongly twisted lines (|T„| > 1.0) are capable of causing X-class flares. Figure [3]^a) shows the 
field lines in orange with the strong current density region represented by the green surface 
over the normal component of the magnetic field in gray scale at 20:30 UT on December 
12. Because these field lines surround the surface of the strong current density region, they 
are expected to store considerable free energy. Figure EJj^b) represents the twist distribution 
by red (|T„| = 0.5) and green (|T„| = 1.0) contours mapped on the photosphere, with the 
same field lines as in Figure ^a). The footpoints of most field lines are clearly rooted in 
the regions surrounded by the red contours, which clarifies that the twist values of the field 
lines supporting the strong current density are distributed in the range of 0.5 < |T„| < 1.0. 
Hence, these results suggest that the field lines having a twist of more than one turn could 
not be deeply connected with the release of free energy in X-class flares. 



3.2.2. An Evolution of the Twisted Lines Near the Polarity Inversion Lines 

In this subsection, we present the temporal evolution of the magnetic twist near the PIL. 
Because many flares are observed to originate at the PIL, it becomes imperative to observe 
the behavior of the magnetic twist near this region. The left panels in Figure H] shows the 
temporal evolution of the twist distribution near the PIL. Blue and red lines indicate the 
locations of the negative and positive sunspots, respectively, which are defined as in Figure 
[U The green lines represent the PIL. The black and white areas represent the negatively 
and positively twisted regions occupied by closed field hues, whereas the regions occupied by 
open field lines are shown in gray and are outside the scope of discussion in this study. At 
17:00 UT on December 11, this active region consisted almost entirely of field lines having 
negative twist values. The positive sunspot's motion had already developed this twisted 
region more than a day before flare onset. The positively twisted regions gradually increased 
as the flare onset approached. Finally, after the flare, most regions near the PIL between 
the positive and negative sunspots were occupied by positively twisted field lines. The time 
evolution of field lines in this AR is characterized mainly by the negatively twisted field due 
to the positive sunspot's motion (Figure H]). However, positive twist begins to buildup at 
03:50 UT on December 12, about one day before flare onset, and it form ed in less time tha n 



the strongly negatively twisted regions. This result is consistent with iPark et al.l ( 120101 ) 



which describes the injection of positive helicity into a pre-existing system having negative 



helicity on a shorter time scale. Recently, iRavindra et al.l ( 120111 ) also showed the injection 
of an opposite vertical current near the PIL through the net current distribution on both 
polarities. 



We also check the average force- free a, denoted as a. The color map in the right panels 
in Figure H] shows the temporal evolution of a. Its pattern is also very similar to that 
of the twisted field lines shown in the left panels. The negative a is already injected in 
the early phase, whereas the positive a increases gradually near the PIL as time proceeds. 
Unlike the twist distribution, some of the positive a, for example, the region surrounded by 
the solid square in the right panels, is comparable in strength to the negative a, although 
the magnitudes of the positive twist values are much smaller than those in the negative a 
regions (see Figure [UJ^b)). Consequently, we see a much different picture in terms of the twist 
analysis. 

3.3. Release Process of the Magnetic Twist in the X3.4 Class Flare Occurred 

in the Solar Active Region 10930 

3.3.1. Twist distribution vs. Evolution of the Flare Ribbon 

Here, we investigate the release process of the magnetic twist before and after the 
flare in detail in order to understand the flare dynamics in terms of the variations in the 
twist value. First, we check the temporal evolution of the Call illumination obtained from 
SOT / Hinode and compare it with the magnetic twist obtained through the NLFFF to clarify 
the relationship between them. 

The gray scale map in Figure [5] shows running difference images of Call from 02:14 
UT to 02:24 UT on December 13. This time period corresponds to the growth phase of 
the two-ribbon structure. The white solid and dashed lines represent the locations of the 
positive and negative sunspots at 20:30 UT on December 12, which are defined as in Figure 
[TJ At the initial onset phase, around 02:14 UT on December 13, the Call image shows 
strong illumination around the local area nearby the PIL, which is surrounded by the white 
dotted square. The two-ribbon structure had still not formed at this time. A thin two-ribbon 
structure started to form around 02:18 UT on December 13 and finally became thick and 
showed strong illumination from 02:22 UT to 02:24 UT on December 13. 

Figure [5] shows contours corresponding to the twist value |T„| = 0.5 (red lines), which 
were obtained from the NLFFF at 20:30 UT on December 12. The regions surrounded by 
the red lines are occupied by strongly twisted field lines (greater than a half-turn twist, 
i.e., \Tn\ > 0.5). The initial illumination at 02:14 UT on December 13 clearly comes from 
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the weakly twisted region where |T„| < 0.25. Although the initial two-ribbon structure at 
02:18 UT on December 13 formed in the weakly twisted region, it seems to be associated 
with strong illumination in the later stage (02:22 UT and 02:24 UT on December 13) in the 
strongly twisted region. From this result, we infer that the strong Call illumination seems 
to be related to the strongly twisted magnetic fields. 

We further investigate the relationship between the Call intensity and magnetic twist 
value more quantitatively. We define the average Call intensity (/) and average twist value 
|(T„)|. The Call intensities are also based on the running difference images, as shown in 
Figure [61 Here "average" indicates the average intensity in the northern ribbon, where only 
the region of strong Call illumination is selected by setting the threshold value of the Call 
intensity to / = 200. However, because the two-ribbon structure had not formed at 02:14 
UT on December 13, we cannot distinguish the northern and southern ribbons. Hence, we 
need to calculate these average intensities using another path instead of the northern ribbon. 
Figure El^a) shows the area surrounded by the dotted square in Figure [5](a); the strong CaH 
intensity (filled contours at / > 150) at 02:14 UT on December 13 is plotted on a distribution 
map of the normal component of the magnetic field (gray scale) at 20:30 UT on December 
12. At this time only, (J) and \{Tn)\ were calculated for the strong Call intensity region 
surrounded by a dotted circle in Figure |6]J^a). Figure |6]^b) shows the profile of (/), which 
is normalized by the (J) = 1000, as a function of |(T„,)|. This result clearly shows that (/) 
suddenly increased for |(T„)| = 0.47 at 02:22 UT on December 13, which indicates that a 
dramatic reconnection occurred in the magnetic twist having values of |(T„,) | > 0.47 because 



strong Call illumination is strongly correlated with magnetic reconnection fjPriest fc Forbes 



2OO2I ). From these results, we find that the strong twist has a significant relationship with 



the strong Call intensity. 



Note that this conclusion is not inconsistent with our previous finding (jinoue et al. 



(120121 )). in which we showed that strong X-ray intensity is not related to strongly twisted 
lines but rather shows a good correlation with the field- aligned current. Our previous study 
examined the core region of the sigmoid, which consisted of weakly twisted field lines and 
was located inside the strongly twisted field lines forming the elbow part of the sigmoid 
and revealing strong Call intensity. Therefore, the environment of the magnetic field lines 
making up the core of the sigmoid differs from the environment we discuss in this work. 



3.3.2. Relaxation of the Magnetic Twist through the Flare 



We proceed with a detailed investigation of the changes in the magnetic topology due 
to magnetic reconnection before and after the flare by using the magnetic twist obtained 
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through the NLFFF and Call images. Figures [7|^a) and (b) show the magnetic field fines (in 
gray) before and after tfie fiare over a time integration of tfie Call images from 02:20 UT 
to 02:40 UT on December 13, tfiat is, from tfie growtfi to tfie formation of tfie two-ribbon 
flare structure. Tfie regions surrounded by red fines, labeled Rl and R2 or Rl' and R2', 
are occupied by closed strongly twisted lines > 0.5) from tfie NLFFF (a) before and 
(b) after tfie flare (20:30 UT on December 12 and 04:30 UT on December 13, respectively). 
Selected magnetic fleld lines, labeled LI (LI') and L2 (L2') connect tfie strongly twisted 
regions Rl (Rl') and R2 (R2') before (after) tfie flare. Note tfiat altfiougfi LI (or L2) and 
LI' (or L2') are plotted for tfie same location, tfiere is no similarity among tfiem. 

Tfie strongly twisted regions Rl and R2 before tfie flare are clearly located in tfie area 
wfiere tfie Call image sfiowed brigfitening, some parts of wfiicfi disappeared after tfie flare 
even tfiougfi regions Rl' and R2' remained as sucfi. In tfie region wfiere tfie strong twist 
disappeared, tfie selected magnetic fluxes LI' and L2' seemed to relax into an untwisted 
fle ld relative t o LI and L2. In particular, tfie fleld line LI forms a compact loop, as sfiown 



m 



Sun et al.l ( 12012| ). Because tfiis relaxation appears on tfie image of tfie strong Call il- 
lumination, magnetic reconnection may be a candidate for explaining tfie relaxation of tfie 
magnetic twist. 

Figures [7](c) and (d) sfiow scatter plots of tfie twist T„ (vertical axis) vs. tfie normal 
component of tfie magnetic fleld (fiorizontal axis) before and after tfie flare, respectively. 
Negatively twisted lines are distributed over a wide range of tfie normal component of tfie 
magnetic fleld. Dasfied circles A and B indicate tfie strong negatively twisted regions in tfie 
negative and positive polarities, respectively, before tfie flare. Many of tfie points in regions 
A and B before tfie flare clearly disappeared after tfie flare, altfiougfi some remained. Tfie 
twisted fleld lines representing T„ ^ —0.5 seem to still be distributed over a wide range of 
tfie normal component of tfie magnetic fleld after tfie flare. Tfierefore, we suggest tfiat most 
of tfie twist release is caused by fleld lines witfi at least < —0.5. On tfie otfier fiand, 
tfie positively twisted regions marked by dasfied circle C appear for \Bz\ < 0.2 before tfie 
flare and remain as sucfi even after tfie flare. From tfiis result, we clearly see again tfiat tfie 
mixed region indicating coexistence of tfie negatively and positively twisted lines is formed 
at \Bz\ < 0.2, as sfiown in tfie left panels in Figure Si 



4. Discussion 

In tfie previous section, we discussed tfie buildup and release processes of tfie magnetic 
twist tfirougfi tfie fiuge X-class flare associated witfi AR 10930. Here, we address a possible 
mecfianism for tfie onset of tfiis flare. 
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The above results (Figures [5] - [7]) showed that energy is released mainly by magnetic 
reconnection of the field lines having a strong negative twist. However, we still have not 
answered the question regarding the fiare onset mechanism that we posed at the beginning 
of the paper. Figure M^a) shows that fiare onset starts at around 02:14 UT between the 
positive and negative polarities of the sunspot, where weakly twisted lines (|T„| < 0.25) 
appear near the PIL. Because this region is also composed of a mixture of field lines having 
both positive and negative twist, magnetic reconnection may be easily induced between lines 
having opposite twist. The mixed-twist region was able to form because positively twisted 
field lines were injected across the pre-existing negatively twisted lines. 

Further, we investigate the temporal evolution of the magnetic fiux corresponding to the 
positive twist and compare it with that of the magnetic fiux corresponding to the negative 
twist. Figure [H] shows the temporal evolution of the magnetic fiux, log ($) (5, > 0), for 
Tn < —0.5 and T„ > 0, which are represented by solid and dashed lines, respectively, 
corresponding to the areas displayed in Figures [2] and HI respectively. The vertical dashed 
line indicates the fiare onset time. We found that most of the magnetic twist producing the 
fiare was greater than a half-turn in the negative sense (T„ < —0.5). The time profiles of 
both magnetic fiuxes show an increase as the fiare onset time approached. However, the 
magnetic fiux related to T„ < —0.5 decreased dramatically around the fiare onset time. We 
believe that this decrement in the magnetic fiux related to the negative twist is due to the 
relaxation via magnetic reconnection shown in Figures [MZl In contrast, the magnetic fiux 
corresponding to T„ > increased continuously before and after fiare onset. 



In previous studies, llnoue et al.l (120111 ) and llnoue et al.l ( 12012| ) indicated that the mag- 
netic configuration occupying the strongly twisted regions before the fiare was robust against 
kink mode instability; therefore, this AR cannot be destabilized through kink instability. 
Hence, fiare onset is triggered via another mechanism: breaking of the equilibrium condi- 
tion of the strongly twisted field; the buildup of positive twist may be a candidate cause 
of that collapse. Consequently, we can say that the increase in the magnetic fiux due to 
injection of twist of the opposite sign across the pre-existing field is an important process 
for understanding the fiare trigger mechanism. 

We again look at the distribution maps of the twisted and average force-free a (denoted 
as a) at 20:30 UT on December 12 in the left and right panels in Figures |H Although 
the twist value in the mixed-sign region marked by the black square in Figure [T]^a) is weak 
relative to the negative twist values distributed at a considerable distance from the PIL 
(Figure lU^b)), the strong a is injected into the region occupied by this weak mixed-sign field 
(for example, the region surrounded by the solid square in the right panels in Figure S]). The 
twist value is small in this region because it depends strongly on the field line length, even 
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though strong force-free a exists at their footpoints. iPevtsov et al.l (Il997( ) reported a strong 
correlation betwee n the force-free a and the geometric shear angle Q (equation (5.3.18) in 



Aschwandenll2004j ) from the soft X-ray loops observed by the Soft X-Ray Telescope (SXT) on 



Yohkoh. Therefore, this result implies that the buildup of strong magnetic shear of opposite 
sign may be related to flare onset, even though the field lines are extremely short (see Figure 
[U^c)). However, we cannot identify the essential process related to the amount of positively 
twisted magnetic flux or the strength of the average force- free a. We need further statistical 
analysis using data with a higher time cadence or numerical experiments to answer this 
question. 

Some authors have already proposed theoretical fiare/CME models in which small- 
scale reconnection in the local area at the lower corona leads to l arge-s cale reconnection, 



releasing the magnetic energy. In a 2D model, IChen fc Shibatal (120001 ) and IShiota et al. 



(120051 ) successfully interpreted an eruption as magnetic reconnection between the emerging 
flux and the pre-existing magnetic flux, which leads to a loss of equilibrium in the flux tube; 
subsequently, another magnetic reconnection is induced in the pre-existing lines under the 
eruptive flux rope. 



On the other hand, iKusano et al.l (120041 ) conducted a 3D simulation to determine the 
role of local instability to determine whether it leads to an ejective eruption. They reported 
that tearing instability occurs at the lower corona where the reversed magnetic shear field is 
formed against the overlying field and eventually leads to ejective eruption through a non- 
linear feedback process, that is the mutual ii iteraction between large-scale strongly twisted 
lines and a small-scale region of mixed twist (IKusano et al.ll2004j ). They also proposed that 
the key process of flare onset is magnetic reconnection between field lines having different 
twist in the local area of the lower corona, which becomes an important cause of large-scale 
fla res. Thus, our result shows behavior similar to that of the theoretical model proposed 
by iKusano et al.l (120041 ). However, it is much difficult to conclude an answer of the onset 
mechanism in term of the only quasi-static pictures obtained from NLFFF. 

Although we suggest that the formation of two different regions is important for flar e 
onset, the detai l ed form ation process is still n o t well known. Nevertheless, iPark et al.l (120101 ) . 
Ravindra et al. (2011), and Ravindra et al. (j2011 ) recently suggested that the process of 



building up positively twisted lines over pre-existi ng negatively twist ed li nes seems to be a 
possible component of the flare onset mechanism. I Wang et al.l (120081 ) and iLim et al.l (120101 ) 
indicated that flux emergence is a key process in the formation of these structures and is also 
responsible for flare onset. The key problems are how the mixed-helicity region is formed on 
the local area near the PIL and how it plays a key role in exciting a flare. A flux emergence 
simulation may address these questions in the future. 
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5. Summary 

In this paper, we investigated the temporal evolution of the 3D magnetic structure of 
AR 10930. We examined the buildup and release processes of the magnetic twist associated 
with this AR using a time series of the vector fields before and after a flare. Our analysis 
is based on the NLFFF extrapolation method, which is applied to time series data obtained 
from SP / Hinode. We also suggested a possible flare onset mechanism in terms of magnetic 
reconnection between field lines having opposite twist. 

We found that the 3D magnetic structure of this AR can be roughly separated into two 
regions containing field lines having different twist. One is occupied by field lines of strong 
negative twist, (z.e., left-handed twist) rooted in the regions of both polarities and is located 
at a considerable distance from the PIL. The other is of mixed polarity and contains both 
positively and negatively twisted lines near the PIL between two sunspots. We investigated 
the temporal evolution of the magnetic twist {i.e., from its buildup to its release) associated 
with this AR using time series NLFFF. We found that a day before flare onset, most regions 
were occupied by negatively twisted field lines. Positively twisted lines were also built up near 
the PIL within one day; consequently, the mixed region was formed, which was eventually 
occupied and dominated by positively twisted lines after the flare. 

For flare onset, we suggest the importance of the buildup near the PIL of twisted 
magnetic lines of different sign compared to the ambient field. In this situation, magnetic 
reconnection can be induced between field lines having different twist, which would destroy 
the equilibrium of the magnetic configuration. We also found that the main relaxation 
process in the flare dynamics is caused by magnetic reconnection in strongly negatively 
twisted lines by comparing Call images with t he NLFFF. This scenario is similar to a 



previous model proposed by iKusano et al.l ( l2004l ). Unfortunately, we cannot extract much 



more information from the NLFFF and conclude a clear answer for the onset mechanism 
because it is useful for magnetic structures in a quasi-static state. Thus, we still have 
problems such as the transition process from a stable condition to an unstable or non- 
equilibrium condition and the 3D dynamics of magnetic reconnection in the solar corona. 
We believe that MHD simulations may provide strong clues regarding these points. 

In addition, this study is essentially limited by the low time cadence of the SP data (~10 
h) used in this analysis. To achieve a more comprehensive view of the temporal evolution 
of the 3D magnetic structure, high-cadence data are essential. A new solar physics satellite. 
Solar Dynamics Observatory (SDO) was launched recently by the National Aeronautics and 
Space Administration (NASA). The Helioseismic and Magnetic Imager (HMI) on board 
SDO can provide vector magnetograms with higher time cadence and a larger field of view 
compared to the data obtained from Hinode. The detailed time evolution of the magnetic 
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field lines and accurate reproduction of the location of the separatrix are important issues 
for the understanding of the flare onset and dynamics. We will deepen our understanding 
of the onset of solar flares by combining high temporal resolution data from SDO with high 
spatial resolution data obtained by Hinode. 
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L[10cm] (c) Bz[kG] (d) 




Fig. 1. — (a) Distribution of the normal component of the magnetic field (Bz) on the 
photosphere at 20:30 UT on December 12 (6 h before the flare). White solid lines indicate 
selected contours of (=790 G and —790 G); blue dotted lines are the PIL. Black solid 
square indicates the region near the PIL between the positive and negative polarities, (b) 
Twist (Tn) distribution of the field lines from NLFFF at 20:30 UT on December 12 mapped 
on the photosphere with the same white and blue contours as in (a). Positive and negative 
signs represent the distributions of right- and left-handed twisted lines, respectively, (c) 
Field line length mapped on the photosphere. Thin white lines have the same meaning as 
in (a); thick white diagonal lines indicate the region occupied by the open field lines that 
are connected outside the field of view. Red lines indicate a twist magnitude of |T„| = 0.5; 
i.e., the regions surrounded by red lines represent |T„| > 0.5. (d) Selected field hues plotted 
over the B^ distribution. Blue, orange, and green field lines correspond to twist ranges of 
\Tn\ <0.5, 0.5< \Tn\ <1.0, and \Tn\ >1.0, respectively. All panels(a— d) have the same field 
of view. 
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Fig. 2. — Upper panels(a— e): temporal evolution of selected twist magnitudes |T„| = 0.5, 
\Tn\ = 1.0, and |T„| = 1.5, plotted in red, green, and blue lines, respectively, over the normal 
component of the magnetic field in gray scale in the same field of view as Figure [1] Lower 
panels(a'— e'): selected magnetic field lines are also plotted. Orange and green field lines 
represent twist ranges of 0.5 < |T„| < 1.0 and 1.0 < |T„| < 1.5, respectively. 
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Bz[kG] (3) Field Lines + Current (b)Field Lines + Twist 




Fig. 3. — Magnetic field lines and surface having strong current density (| J| > 2.4) at 20:30 
UT on December 12 (orange and green, respectively) over the normal component of the 
magnetic field (gray scale) in the same field of view as Figure [H (b) Twist distributions in 
red {\Tn\ = 0.5) and green (|T„| = 1.0) contours with the field lines as in (a). 
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Fig. 4. — Left panels; Positive and negative twist distributions for closed loops plotted 
in white and black, respectively. Gray regions are occupied by open field lines. Red and 
blue lines indicate the locations of the positive and negative sunspots, respectively, which 
correspond to a normal component of the magnetic field of 790 G and —790 G, respectively. 
Dotted green lines indicate the PIL. Right panels; Average force-free a, denoted by a, 
plotted in color. White solid lines indicate positive and negative polarities; white dotted line 
represents the PIL, which is also shown in (a). 
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I (a) 02:14 UT- 02:12 UT (b) 02:18 UT - 02:16 UT 




(c) 02:22 UT - 02:20 UT (d) 02:24 UT - 02:22 UT 




Fig. 5. — (a)-(d) Gray scale maps showing the temporal evolution of Call images, which are 
running difference images, corresponding to growth of the flare ribbon in the early phase of 
the flare (02:14 UT to 02:24 UT on December 13). The field of view is the same as in Figure 
[H Red lines indicate contours of twist |T„| = 0.5 obtained from the NLFFF at 20:30 UT 
on December 12. White solid and dashed lines represent contours of the normal component 
of the magnetic field, defined as in Figure dJ^a). Dotted square suggests the location of 
illumination corresponding to the flare onset. 
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Bz[KG] (a) 02:1 4 UT - 02:1 2 UT (b) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 
|<Tn>| 

Fig. 6. — (a) Areas of strong Call intensity (/ >150, filled contours) at 02:14 UT on 
December 13 plotted on a distribution map of the normal component of the magnetic field 
(gray scale) at 20:30 UT on December 12. Image covers the area enclosed by the dotted 
square in Figure [5](a). Northern and southern ribbons cannot be distinguished at this time 
because the two-ribbon structure had not formed yet, as shown in Figure O (b) Profile of 
(/), normalized by the (/) = 1000, as a function of at various times on December 13. 
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Fig. 7. — Color maps of time integration of Call images from 02:20 UT to 02:40 UT on 
December 13. Red lines represent twist values of |T„| = 0.5 from the NLFFF at (a) 20:30 
UT on December 12 (before the flare) and (b) 04:30 UT on December 13 (after the flare). 
Rl (Rl') and R2 (R2') represent the strongly twisted regions {\Tn\ >0.5) on the negative 
and positive polarities, respectively. LI (LI') and L2 (L2') are selected magnetic field lines 
connecting Rl (Rl') and R2 (R2'). (c) and (d) Scatter plots of twist T„ (vertical axis) vs. Bz 
(horizontal axis) before and after the flare, respectively. Dashed circles A and B in (c) show 
the regions where the strong twist was distributed before the flare. Region C represents the 
positively twisted region of approximately \Bz\ < 0.2 before the flare. 
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Fig. 8. — Temporal evolution of the magnetic flux, log(<l') {B^ > 0), for T„ > (dashed 
line) and T„ < —0.5 (solid line), respectively, corresponding to the areas displayed in Figures 
|2] and HI respectively. These values are normalized by 10^^ (Mx) (= J B^dS; dS represents 
an surface element in the bottom surface). Vertical dashed line indicates the time at which 
the flare occurred. 



